Broken-Symmetry States of Dirac Fermions in Graphene with A Partially Filled High 

Landau Level 
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We report on numerical study of the Dirac fermions in partially filled N = 3 Landau level (LL) in 
graphene. At half-filling, the equal-time density-density correlation function displays sharp peaks 
at nonzero wavevectors ±q*. Finite-size scaling shows that the peak value grows with electron 
number and diverges in the thermodynamic limit, which suggests an instability toward a charge 
density wave. A symmetry broken stripe phase is formed at large system size limit, which is robust 
against purturbation from disorder scattering. Such a quantum phase is experimentally observable 
through transport measurements. Associated with the special wavefunctions of the Dirac LL, both 
stripe and bubble phases become possible candidates for the ground state of the Dirac fermions in 
graphene with lower filling factors in the N = 3 LL. 

PACS numbers: 73.43.-f, 71.10.-w, 73.22.Gk, 71.45.Lr 



Recently, successful fabrication of single-atomic-layer- 
thick films of graphite [![, called graphene, has triggered 
intensive research activities to understand the novel prop- 
erties @, H, S S Hi of these new two-dimensional (2D) 
electron systems (2DES's). Different from the conven- 
tional 2D electrons with the quadratic dispersion rela- 
tion, the low-energy electron excitations in graphene have 
a linear (relativistic) dispersion relation, which can be 
described by a massless Dirac equation @, HI, 0- The 
Dirac-fermion-like nature of the electrons has manifested 
itself evidently in the unconventional quantization pat- 
tern of integerquantum Hall effect (IQHE) 0,0, 0, |E 
HI, EH EH, EH E3 • The possibilities to observe the frac- 
tional quantum Hall effect in graphene in the partially 
filled N = and N = 1 LLs and a pseudospin ferro- 
magnetic state in the N = 1 LL have been predicted in 
some recent theoretical works [IH El El El, which how- 
ever have not been experimentally observed yet. Thus 
the role of impurity scattering seems crucially impor- 
tant, which may obscure the observation of these quan- 
tum phases for currently available sample mobility. At 
even higher LLs, there exist many studies on the con- 
ventional 2DES's, where a variety of quantum phases 
have been predicted |17I . Il8l . |19H and discovered ex- 
perimentally [20( | . varing from stripe phase to reentrant 
QHE. On the other hand, for the Dirac fermions with 
partially filled high LLs in graphene, investigations are 
only carried out based on the Hartree-Fock mean-field 
approximation [2l| , which suggested that stripe and bub- 
ble phases are possible. It remains an open issue whether 
such quantum phases can survive quantum fluctuations 
and disorder scattering. Study of the quantum phases of 
Dirac fermions in high LLs is highly valuable for eluci- 
dating the key role of the Coulomb interaction and the 
disorder effect in graphene. 

In this Letter, we investigate the low-energy states 



of the Dirac fermions in graphene with partially-filled 
A^=3 Dirac LL using Lanczos method for finite-size sys- 
tems with torus geometry and upto N e = 16 electrons 
in the LL 22]. For a pure system at half- filling, we find 
a large number of nearly degenerate low-energy states, 
which are equally distanced from each other by a charac- 
teristic wavevector q* in momentum space. The equal- 
time density-density correlation function shows strong 
and sharp peaks at ±q*, indicating the charge density 
wave (CDW) instability. Finite-size scaling shows that 
the peak value grows with electron number and becomes 
divergent in the thermodynamic limit. Thus a symmetry 
broken stripe phase will be formed at large system size 
limit, which is robust against moderate disorder scat- 
tering with a critical disorder strength comparable with 
that for the nonrelativistic 2DES's [23J]. Due to the spe- 
cial single-particle wavefunctions of the Dirac LL, both 
stripe and bubble phases are possible candidates for the 
ground state with lower filling factors in the N = 3 LL 
based on exact calculations. We further discuss the trans- 
port anisotropy of the stripe phase, which can be used to 
experimentally detect this quantum phase. 



We consider a 2DES in an L x x L y rectangular cell 
of graphene under a perpendicular magnetic field. The 
magnetic length I is taken to be the unit of length. The 
total number of flux quanta = L x L v /2tt is chosen to 
be an integer. Periodic boundary conditions are imposed 
in both x and y directions [l8j]. The magnetic field is as- 
sumed to be strong enough so that the spin degeneracy of 
the LLs is lifted, and all the LLs are well separated from 
each other. One can thus project the system Hamilto- 
nian into the topmost, partially filled, N-th LL [18]. The 
projected Hamiltonian, which contains the Coulomb in- 
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FIG. 1: Low-energy spectrum for the half-filled N = 3 LL 
with N e = 12 in a clean system, (a) Energy levels versus 
aspect ratio asp, where the eigenenergies are measured in 
units of e 2 /e£ with respect to the ground state, (b) Energy 
levels versus the y component of wavevector at aspect ratio 
asp = 0.74. The six nearly degenerate lowest-energy states 
are equally separated by a wavevector q* = (0,4n/L y ) — 
(0,0.88). 

teraction and disorder potential, has the form: 

h c = j2J2 e ~ q2/2 ^^ 2v ^y MRi ' Rj) / 27rN ^ 

i<j q 

+ EE eV/4 ^(9)W<zK q - Ri , (i) 

i q 

where is the guiding center coordinate (GCC) of the 
i-th electron, and V(q) = 2ire 2 /eq is the Fourier trans- 
form of the Coulomb interaction. The wavevector q takes 
discrete values that are compatible with the geometry of 
the system. Noticing that the electron wavefunction in a 
relativistic Dirac LL of N is a mixture of two wavefunc- 
tions in two different nonrelativistic LLs of N and N — 1, 
we can write the form factor Fjv(q) in Eq. (1) as [13] 

F N (q) = l -[L N {q 2 /2) + L N _ 1 {q 2 /2)] , (2) 

where Ln(x) is the Laguerre polynomial. The disorder 
potential is generated according to the correlation rela- 
tion in g-space (V imp (q)V imp (-q')) = j^Sq-q', which 
corresponds to (Vi mp (r)V r i mp (r')) = W 2 S(r — r') in real 
space with W as the strength of disorder in units of e 2 /e£. 

We compute exactly the low-energy spectrum and 
wavefunctions using the Lanczos diagonalization method. 
We find that pseudospin-polarized states have lower en- 
ergies than unpolarized states. In Fig. la, the low-energy 
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FIG. 2: (Color online) Correlation functions of the ground 
state for the half-filled N = 3 LL in a clean system, (a) Static 
density-density correlation function •S'o(q) with N e = 12 and 
asp — 0.74. (b) Peak value of So(q) normalized by N e as a 
function of 1/N e for iV=2 and 3 LLs. (c) Ground-state pair 
correlation function g(r) in guiding center coordinates with 
N e = 12 and asp = 0.74. 



spectrum for the N — 3 LL with filling factor v — 1/2 
and N e = 12 is shown as a function of the aspect ra- 
tio asp — L x /L y . Here, the filling factor is defined as 
v = N e /N<f, with N e the electron number in the LL. A 
generic feature of the spectrum is the existence of a num- 
ber of nearly degenerate low-energy states well separated 
from higher-energy states by a gap, in a wide range of as- 
pect ratio 0.7 < asp < 0.9. We will call these states the 
ground-state manifold. In Fig. lb, we plot the low-energy 
spectrum as a function of wavevector q y at the optimized 
aspect ratio asp — 0.74, where the energy broadening of 
the manifold is minimized. All the states in the manifold 
are equally separated from each other by a characteristic 
wavevector q* = (0,q*) with q* — 4n/L y . This feature 
is similar to that observed in the spectra of the non- 
relativistic 2DES's with half-filled high LLs [3, which 
indicates an instability towards a unidirectional CDW or 
stripe phase. 

We next turn to the LL projected equal-time 
density-density correlation function S'o(q) of the ground 
state |li]. In Fig. 2a, 5 (q) for the half-filled N = 3 LL 
with N e = 12 and aspect ratio asp = 0.74 is shown in 
a 3D plot. We see that there are two sharp and strong 
peaks at ±q* = (0, ±0.88) with the peak value 4.37, while 
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FIG. 3: Disorder effect on the ground state of the half-filled 
N = 3 LL with iV e = 12 and asp = 0.74. Projected elec- 
tron density p(r) at disorder strengths (a) W = 0.01, (b) 
W = 0.10, and (c) W = 0.20. (d) Peak value of 5 (q) versus 
disorder strength. The transition occurs around W = 0.16. 



So(q) is about 0.5 away from the peaks. The presence of 
the peaks in Sb(q) at q* suggests strong density correla- 
tion at the ordering wavevector, which is consistent with 
the characteristic feature of the low-energy spectrum. We 
have also examined So (q) at a number of other aspect ra- 
tios in the range of 0.7 < asp < 0.9, which all show sharp 
peaks but with slightly reduced peak values. The ratio 
So j^ - for different electron numbers from N e = 6 upto 
N e = 16 is shown as a function of 1/N e in Fig. 2b, at the 
corresponding optimized aspect ratios. The ratio 5 °^ q ^ 
extrapolates linearly to a finite value 0.28 as iV e — > oo, 
which is proportional to the relative density modulation 
in the symmetry broken state. Therefore, the instability 
towards a unidirectional CDW phase in large systems is 
established by finite-size scaling. We have also studied 
the iV=2 LL, where both the value of So(q*) and the 
value of at N e — > oo are smaller than those in the 

N = 3 LL case suggesting the N = 3 LL to be the better 
candidate for observing the CDW phases in graphene. 

The CDW order of the above ground state in the real 
space is studied using the LL projected pair correlation 
function g(r) 



24j . In Fig. 2c, we plot the ground- 



state pair correlation function in GCC for the half-filled 
N = 3 LL with N e = 12 and asp = 0.74. It shows 
clearly that there are two stripes inside the unit cell along 



the x direction. The mean separation between the two 
stripes is related to the ordering wavevector q* through 
D s = 2ir/q* = 7.1, which happens to be a half of L y 
for the present system. By examining systems with N e 
ranging from 6 upto 16 at their optimized aspect ratios, 
we find that the ordering wavevector q* changes slightly 
with N e between q* — 0.88 and 0.96, and correspond- 
ingly D s varies between 6.5 and 7.1 magnetic lengths, 
which are comparable to the ones for the nonrelativistic 
2DES's 
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The above calculated correlation functions establish a 
stripe phase at the thermodynamic limit, in agreement 
with the Hartree-Fock result [2lJ . We also note that in a 
clean finite-size system (W = 0), the ground state wave- 
function is invariant under the magnetic translation of 
guiding centers of the cyclotron orbits, while the sym- 
metry broken stripe phase can become the true ground 
state when a relatively weak disorder is turned on [23j |. 
In the presence of such random disorder, we can directly 
probe the CDW phase using the LL projected local elec- 
tron density p(r) For a very weak disorder strength, 
the typical behavior of p(r) is shown in Fig. 3a, where 
W = 0.01. Two nearly perfect stripes are formed along 
the x direction with essentially no density modulations 
along the stripes. At a intermediate disorder strength 
W = 0.1, there are still two complete stripes, though 
the shape of the stripes exhibits pronounced variations 
along the stripe direction. Thus we conclude that for 
weak to moderate disorder strength, the ground state is 
in the quantum Hall stripe phase. With further increas- 
ing disorder strength, as shown in Fig. 3c for W = 0.2, 
the stripe density profile becomes raptured and riddled 
with defects, indicating that the long-range stripe order 
is lost. The remaining short-range stripes become ran- 
domly orientated. In the relatively weak disorder region 
< W < 0.1, the whole structure of the density-density 
correlation function So(q) is found to be similar to the 
clean system case, indicating the robustness of the CDW 
order. In Fig. 3d, we plot the peak value So(q*) as a 
function of W . The peak value So(q*) remains nearly 
constant for relatively weak disorder W < 0.1. It starts 
to drop quickly around W — 0.13 and becomes compa- 
rable with the background value 0.5 at W > 0.16. Thus 
W = 0.16 is determined as the critical disorder strength, 
where the long-range CDW order becomes unstable. 

The transport property of the above anisotropic CDW 
state is studied by calculating the Thouless energy AE TT 
(r = x or y), i.e., the ground state energy difference due 
to a change of the boundary condition in the r direction 
from periodic to antiperiodic, which is related to the lon- 
gitudinal conductance in the x or y direction. We find 
that for a weak disorder, e.g., W — 0.02, AE yy is about 
50 times smaller than AE XX , suggesting a large trans- 
port anisotropy associated with the quantum Hall stripe 
phase, which can be used to identify the stripe phase 
exp er iment ally. 
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FIG. 4: Ground-state pair correlation functions for the N = 3 
LL for filling factor v = 1/3 and electron number iV e = 8 at 
the aspect ratios (a) asp = 0.65 and (b) asp = 0.86. 

We further study the ground state at lower filling fac- 
tors in the N — 3 LL. For v — 1/3 and N e = 8, the low- 
energy states are found to become nearly degenerate in 
two different aspect ratio regions. The ground-state pair 
correlation functions corresponding to these two regions 
are plotted in Fig. 4(a) and Fig. 4(b), where the aspect 
ratio is taken to be asp — 0.65 and asp — 0.86, respec- 
tively. In the region around asp = 0.65, the ground state 
is a bubble- like CDW phase [19|] . The bubbles distribute 
in a superlattice structure with the lattice vectors Al and 
A2. For asp = 0.65, we have Al = A2 = 6.24, and the 
angle between the two lattice vectors is 9 = 75°. Within 
the unit cell, there exist four bubbles and each bubble 
contains two electrons. In the region around asp = 0.86, 
the ground state is a stripe phase, and the stripe sepa- 
ration is about D s — 6.62 for asp — 0.86. These results 
indicate that both the two-electron bubble phase and the 
stripe phase can be separately realized, depending on the 
geometry of the system. This is different from the case 
for nonrclativistic 2DES's, where usually only a single 
CDW phase exists at a given filling factor in a given 
LL [l9| • More competing phases in the Dirac LL may be 
explained as a result of the special structure of the Dirac 
LL, which is essentially a combination of the N = 2 and 
N = 3 nonrelativistic LLs. 

In summary, we have shown that various CDW phases 
can be realized in the partially-filled N = 3 LL in 
graphene, the character of which can be explained from 
the mixture property of the electron wavefunctions of the 
Dirac LL. The CDW phases are robust against moder- 
ate disorder scattering. The unidirectional CDW stripe 
phase at half-filling shows large transport anisotropy, 
which can be measured experimentally. 
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